Breakdown of the Fluctuation-Dissipation Theorem for fast superdiffusion 
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Abstract 

We study anomalous diffusion for one-dimensional systems described by a generalized Langevin equation. We show 
that superdiffusion can be classified in slow superdiffusion and fast superdiffusion. For fast superdiffusion we prove that 
the Fluctuation-Dissipation Theorem does not hold. We show as well that the asymptotic behavior of the response 
function is a stretched exponential for anomalous diffusion and an exponential only for normal diffusion. 
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Since its formulation, the Fluctuation-Dissipation 
Theorem (FDT) has played a central role|Q, in non- 
equilibrium statistical mechanics (NESM). It reaches 
such an importance that a full formulation of NESM 
is given Q based on it. In the last 30 years, funda- 
mental concepts and methods have been developed 
and a large number of connections have been 
established (see ref. and references therein). A 
necessary requirement for the FDT is that the time- 
dependent dynamical variables are well defined at 
equilibrium. The presence of far from equilibrium dy- 
namics may lead to situations where the FDT does 
not hold, the aging process in spin-glass systems be- 
ing a good example 

Diffusion is one of the simplest processes by which a 
system reaches equilibrium. For normal diffusion, the 
process is so well known that it may be described by 
an equilibrium type distribution for the velocity and 
position of a particle. However, the strange kinet- 
ics of anomalous diffusion, intensively investigated in 
the last years ||-||l3|], shows surprising results. Con- 
sequently, studying anomalous diffusion seems to be 
the best way to obtain the conditions of validity for 
the FDT. 

In this letter, we present a straightforward proof 
of the inconsistency of the FDT for a certain class 
of superdiffusive processes described by a generalized 
Langevin equation (GLE). The use of the FDT allows 
us to classify two classes of superdiffusion. The first 
class, which we shall call slow superdiffusion, does 
obey the FDT; the second class, which we shall call 
fast superdiffusion, does not obey the FDT. The proof 
is simple and we discuss as well how the diffusive 
process leads to an equilibrium. 

We shall start writing the GLE for an operator A 
in the form [g, §, (§ 



dA(t) 
dt 



T(t-t')A{t')dt' + F{t), (1) 



where F{t) is a stochastic noise subject to the condi- 
tions (F(t)) = 0, (F(t)A(0)) = and 



C F {t) =< F(t)F(0) >=< A 2 > eq T{t). (2) 



Here Cp(t) is the correlation function for F(t) and 
the brackets <> indicate thermal average. Eq. (||) 
is the famous Kubo FDT and it is quite general. In 
principle, the presence of the kernel T(t) allows us to 
study a large number of correlated processes. 

We may naively expect that, by Eq. (|l|) and Eq. 
(□), a system will be driven to an equilibrium , i.e. 



lim < A 2 {t) >=< A 2 > 



eq 



(3) 



We shall see however that this is not always the case 
for superdiffusive dynamics. Let us define the vari- 
able 



y(t) 



with asymptotic behavior 



A(f)dt', 



lim < y 2 {t) >~ t^. 



(4) 



(5) 



For normal diffusion /i = 1, we have subdiffusion for 
fi < 1 and superdiffusion for [i > 1. Notice that if 
A(t) is the momentum of a particle with unit mass, 
y(t) is its position. Using Kubo's definition of the 
diffusion constant we get p3| 



D 



,. < A 2 > eq 
lim — = -. 

*-»<> r(z) ' 



(6) 



where T(z) is the Laplace transform of T(t). A finite 
value of r(0) corresponds to normal diffusion, T(0) = 
to superdiffusion and T(0) = oo to subdiffusion. 
Notice that 



7 = r(o) = / T(t)dt 
Jo 



(7) 



plays the same role as the friction in the usual 
Langevin's equation, i.e., GLE without memory. 
Now we propose a solution for Eq. ml) as 



1 



A(t) = / R(t- t')F{t')dt 
Jo 

where we have set ^4(0) = and (12) 



R{z) 



z + T(z) 



(8) 



(9) 



Squaring Eq. (||) and taking thermal average we ob- 
tain 



Figure 1: Normalized mean square velocity as a function 
of time for the memory given by Eq. (|l9|) . Here (3 — W2/2 
and W2 = 0.5. Each curve corresponds to a different value 
of w\. a) ioi = 0; b) 101 = 0.25; c) Wi = 0.45. The 
horizontal lines correspond to the final average value A s . 
In agreement with the theoretical prediction, A s decreases 
as wi grows. 



b ^ 00 and does violate the FDT. The first class we 
shall call slow superdiffusion (SSD) and the second 
class fast superdiffusion (FSD). 

Consider now the asymptotic behavior for T(z) 



<A 2 (t)>= / C F {t' -t")R(t')R(t")dt'dt". 
Jo Jo 

(10) 

At this point, it is quite usual to perform numerical 
calculation p2| . From Eq. (|l|), we can get a self- 
consistent equation for R(t) as 



= - / T{t-t')R(t')dt'. (11) 
Jo 



dR(t) 
dt 



By using the FDT Eq.(f2|) and Eq.([llJ) we can exactly 
integrate Eq. (10) and obtain 



< A 2 (t) >=< A 2 > eq X(t), 



(12) 



lim T(z) = az 1, 

z^0 



(18) 



For v < 1 we have subdiffusion, for v = 1 normal 
diffusion. For 1 < v < 2 the process belongs to the 
SSD and, finally, for v >2 we have FSD. There is an 
obvious connection between v and /x. Using Eq. (||) 
and the fact that Xmx z ^QY{z) = lim^oo r(l/t) we 
get v — /i and consequently the FSD starts at \i > 2, 
i.e., the ballistic motion and beyond. It is interesting 
to note that Lee || proved the failure of ergodicity 
for the ballistic motion and now we showed that the 
FDT does not hold for this motion. 

Now we test our analysis against simulations. Let 
us consider the function 



where 

A(t) = 1-R 2 (t). (13) 
Notice now that Eq. (H) is satisfied if and only if 



lim X(t) = A* = 1, 

t — >OG 



or equivalently 



lim R(t) = lim zR(z) = 0. 



(14) 



(15) 



Equation ( |l5|) is the ergodic condition ||, It is sat- 
isfied for normal diffusion and subdiffusion. Now for 
superdiffusive systems 



lim R(t) = (1 + b)-\ 

t—*oo 



(16) 



where 



b = lim 



df(z) 



dz 



There are two distinct limits for b, which define two 
classes of superdiffusion. For the first class, b = oc 
and the system obeys the FDT. The second class has 



r(t) = (5 



sin(w2t) sin(wii) 



t 



t 



(19) 



where w 2 > W%. This function was chosen so that 
f(0) = for any w± ^ . Thus, for W\ — we 
have normal diffusion and for any wi / we have 
superdiffusion with [i — 2. If we let (3 — W2/2 we get 
A* as 



A* = 1 — 



2wi 



tl'l + w 2 



(20) 



Any value of A* different from 1 shows the inconsis- 
tency of the FDT in Eq. (||), because we start sup- 
posing the existence of an equilibrium value < A 2 > eq 
and, after an infinite time, we end up with < A 2 > 



A*. No matter the < A 2 > eq that we input in Eq. (|2j), 
we never reach it, except for the trivial null value. 

Now we select A(t) = v(t), the particle's velocity, 
so that < v 2 (t) >—< v 2 > eq X(t). We simulate the 
GLE for a set of 10,000 particles starting at rest at 
the origin and using the memory in Eq. ( |l9| ) with 
W2 — 0.5 and different values of W\, The results of 
these simulations are shown in Fig. 1, where we plot 
(17) < v 2 (t) >. We used the normalization < v 2 > eq = 1, 
so that < v 2 (t) >= A(i), Notice that A(i) does not 
reach a stationary value, rather it oscillates around 
a final average value X s . This value of X s should be 
compared with A* obtained from Eq. 
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Figure 2: A* as a function of the parameter Wi. Each dot 
corresponds to a value of A a obtained from simulations 
like those described in Fig. 1. The line corresponds to 
the theoretical prediction given by Eq.(Eol). 



In Fig. 2 we plot A* as a function of w\ as in 
Eq. (p0|) with W2 — 0.5. We also plot the final average 
values A s obtained from simulations for different val- 
ues of W\. Note that simulations agree with theory 
and X s — > 1 when wi — > . 

For \x > 2, the FSD cannot be described by the 
methods we discussed here. Once the FDT does not 
work, the GLE and the FDT together predict results 
such as null dispersion for the dynamical variable, 
i.e. < A 2 (t — > oo ) >= 0. Moreover, the exponent 
/Lt can be put as fi = 2/ Dp, where Dp is the fractal 
dimension jTJj. Consequently fi > 2 leads to Dp < 1, 
which is not a full curve, but a set of points such 
as the Cantor set, and cannot represent a classical 
trajectory. 

At first sight, the results presented here seem 
strange. Why does the FDT not work for the FSD? 
As we remarked before, 7 in Eq. (Q) plays the same 
role as the usual friction in the Langevin Equation 
that yields R(t) ~ exp(— jt) with a relaxation time 
t = 7 _1 for large times. For both SSD and FSD, 
7 = 00 and the system should not reach an equi- 
librium. 

Now we address the previous question in another 
way: "Why does the FDT work for the SSD? Is it 
really r = r(0) _1 the relaxation time?". In order to 
answer this question one needs to know the asvmp- 
totic behavior of R(t) as t — > 00. From Eq. ([□]) we 
may write 



In R(t) = -T(t) R(t')dt' -tf(z). (21) 
Jo 

In the limit when t ^00 or, equivalently, z = 1/t — > 
0, it is possible to eliminate the first term at the right 
of Eq. (Ell) by using 



I = lim r(t) / R{t')dt' = lim 

t—>oc Jq 



*r(*) 

z + f(z) 



(22) 



Notice that for f (z) = az^ 1 and fi > 0, I — * 
and we get the asymptotic behavior 



lnR(t) = -t r(t')dt' = -tr(0). (23) 
Jo 



The limit in Eq. (|23|) is quite clear for normal dif- 
fusion, where 7 = f (0) is finite, and for subdiffusion, 
where F(0) — ► 00. However, for superdiffusion one 
must look carefully since f (0) — > 0. We use T(z) as 
in Eq (TL8) to obtain 



lim tf(0) = if (1/t) = at 2 



(24) 



We see that Eqs. (g3J) and (|2J) yield R(t -> 00) = 
only for \l < 2, what includes the subdiffusion, the 
normal diffusion and the SSD. For the FSD, \i > 2 
and we shall use Eq. ( |l6| ) to obtain the infinite limit. 
Thus, in this limit process, there is an infinite relax- 
ation time r = 7 _1 for superdiffusion. However, this 
relaxation time can be seen only as a result of an 
evolution, which, for the SSD, is never of the same 
order of t in the limit t — ► 00. Consequently, for long 
times, the SSD presents a finite relaxation time. In 
short, the SSD has in common with normal diffusion 
and subdiffusion the fact that they have a finite re- 
laxation time and obey the FDT. 

Now we can look beyond the exponential aspect of 
the asymptotic solution Eq. ( p3| ) and use Eq. ( p^ ) to 
obtain 



where 



lim R(t) — exp 

t^QO 



= 2- p. 



(25) 



(26) 



For fi ^ 1, r = a^ 1 ^ and for (i = 1, r = 7 _1 = 
f(0) _1 . The function Eq. (|25|) is a stretched expo- 
nential and we shall discuss that in detail below. 

We have important results. First, we obtain a 
stretched exponential associated with anomalous dif- 
fusion, i. e. both subdiffusion and SSD. Also, we 
obtain the exponent directly, not by fitting nor 
simulations, with no reference to a specific system. 
Finally, we show that the relaxation time of the cor- 
relation function is r(0) _1 only for normal diffusion. 
For that case, the correlation function decays as an 
exponential. For subdiffusion and for SSD the re- 
laxation time is associated with the coefficient of the 
main term of f (z) in the limit when z — > 0. Thus 
we can define a relaxation time for both normal and 
anomalous diffusion in the form 



t = lim 



^-"t{z) 



(27) 



Notice that for // = = 1, r = r(0) _1 as expected 
for normal diffusion. 

Let us discuss the very particular behavior of \l — 
0, i.e. the "no diffusion at all" behavior. This can be 
easily obtained by the constant memory r(t) 
which yields for the friction force in Eq. (Q) — 
This is precisely an harmonic oscillator, which does 
not dissipate nor diffuse at all. For this system, we 
have T(z) = ujqZ^ 1 , and R(t) can be exactly solved 
as a cos(ujot) type behavior. As expected, R(t) has 



— "0, 
muj 2 y. 
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no relaxation time. However, using f(z) on Eq. (pTj) , 
we get t = ojq 1 , which is the time scale of the oscilla- 
tion, i.e. the inverse of the frequency. Consequently, 
even in an extreme situation where we do not have a 
relaxation time, Eq. ( p7| ) yields the right time scale 
of the system. 

The research on the striking universality proper- 
ties of slow relaxation dynamics in glass Jg, Q, su- 
percooled liquids pi , liquid crystal polymer jjiql and 
disordered vortex lattice in superconductors p7| has 
been driving great efforts in the last decades. A large 
and growing literature can been found where the non- 
exponential behavior (stretched exponentials) has 
been observed in correlation functions [ |l5[ |l7| . Those 
have in common the fact that they are subject to 
an anomalous diffusion. Peyrard [Q made a model 
for two-dimensional water and, by using Monte Carlo 
simulation, obtained the correlation function with an 
exponent 0.3 < j3 < 0.6. When the temperature de- 
creases, he suggests that (3 — > 1. Using his data in 
Eq. (|26|) , we get (3 ~ 0.75. It would be too naive to ex- 
pect that our simple unidimensional, linear approach 
would describe all the range of complex structures. 
Nevertheless, it may bring an insight to guide us in 
such situations. 

In conclusion, we discussed the stationary behav- 
ior for the mean square value of a dynamical vari- 
able A(t) and noticed that the superdiffusive motion 
must be classified in slow superdiffusive (SSD) and 



fast superdiffusive (FSD). The FSD motion shows an 
inconsistency between the GLE and the FDT. The 
FSD has infinite relaxation time, and consequently 
never reaches equilibrium. This kind of superdiffu- 
sion in which < A 2 (t) >^ with /i > 2 is com- 
mon in hydrodynamical processes. It is not surprising 
that these processes will be far from equilibrium and 
violate the FDT. We pointed out here how it hap- 
pens and precisely where the FDT breaks down. As 
we have already mentioned, spin glasses seem to be 
a rich field for studying these phenomena. Indeed 
experimental || and theoretical works |(| 0] have 
been reported in this area, confirming the violation of 
the FDT. As well, the stretched exponential behav- 
ior found in noncrystaline material is connected here 
with anomalous diffusion. It would be very helpful 
if the exponent \i for those diffusive processes could 
be measured. Another related phenomenon is the 
anomalous reaction rate, which we expect to discuss 
soon. Although anomalous diffusion remains as a sur- 
prising phenomena, we hope that this work will help 
in the centennial effort to understand diffusion and 
the relation between fluctuation and dissipation. A 
generalization of the FDT to include the FSD is nec- 
essary, what will require a deeper understanding of 
systems far from equilibrium. 
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